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The properties of V838 Mon in 2002 November 
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Abstract. We present the results of modelling the 0.45-1 /im spectral energy distribution of V838 Mon for 2002 
November. Synthetic spectra were calculated using the NextGen model atmospheres of Hauschildt et al. (1999), 
which incorporate line lists for H2O, TiO, CrH, FeH, CO, and MgH, as well as the VALD atomic line list. Fits 
to the observed spectra show that, in 2002 November, the effective temperature of V838 Mon was approximately 
2000 ± 100 K. Our theoretical spectra show a comparatively weak dependence on logy. Preliminary analysis of 
the hot star observed together with V838 Mon shows it to be a normal B3V dwarf. 
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1. Introduction 

V838 Mon (also known as EQJ0704.0-0350, GSC 04822- 
00039, IRAS 07015-0346, USNO-A2.0 0825-03833116, 
AAVSO 0659-03) has been the subject of intense in- 
terest sinceits discovery as an eruptive variable by 
iBrown I l|2002|) on January 6, 2002. Two further max- 
ima on t he light curve subsequently occurred, in 2002 
February llMunari et al. 11200^: iKimeswenger et af~l l2002: 
ICrause et al. 1120031) and the optical flux in the V band 
increased by ~ 9 mag. The luminosity peaked at V ~ 7 
in 2002 February. Later, a gradual decline in the V flux 
be gan which, by 2003 Ja nuary, had faded by 8 mag. 

Henden et al~l {2002) reported the presence of a light 



echo subsequently s tudied in detail bv lBond et al. 1 1)2003 
and Cra use et al. 1 ll2005lh The re flecting dust may be 
interstellar llTvlenda et all l2005h or may have origi- 
nated hi_the_envek2pe of_V8 3 8 Mon and ejected in the 
past (|van Loon et al. lEoO^l . Van Loon et al. (2004) re- 
ported the discovery of the multiple dust shells around 
V838 Mon. Alternatively, these shells arose through heavy 
mass loss during the evolution of a massive progenitor 
(Munari et al. 2005). 

The distance to V838 Mon c an be estimated from 
the evolution of the light echo ijHenden et al. I l2002(h 
According to works based on Hubble Space Telescope 



data, its distanc e is > 6 kpc l|Bond et al. 1l2003h . 8±2kpc 
l|Tvlendall2004|) . or as high as 12 kpc (Crause et al. 2005). 
However, according to recent work based on HST data, 
its distance is 5.9 kpc (Sparks et al. 2006). If these esti- 
mates are correct, then at the time of maximum brightness 
V838 Mon was the most luminous star in the Galaxy. 

The cause of the eruption of V838 Mon and the na- 
ture of its progenitor are unclear. The integrated colors of 
the progenitor looke d that of an unredden e d F-t ype star 
l|Munari et al. Il2002|) . iDesidera fe Munari I l)2002l) discov- 
ered spectroscopically a hot companion to the outbursting 
star, later confirmed by Wagner et al. (200 2), and classi- 
fied as a B3V star bv lMunari et al. I l)2005|h 

The spectral, photometric, and polarimetric evo- 
lution of V838 Mon ha s been d escribed in several 
studies 
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of TiO had appeared in the spectrum and by May, 
the spectrum had evolv ed to a "very cold" M giant 
l|Baneriee fe Ashok lEoO^l . Its infrared spectrum in 2002 
October was char acterized as an "L-supergiant" by 
lEvans et all l|2003|) . 

The complexity and the rarity of V838 Mon opti- 
cal spectrum resulted in only a few spor adic attempts 
to der ive the star's atmospheric properties. iKipper et al. I 
(200-iJ) found that, for the iron group of elements, 
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[M/H] = —0.4, while abundances of lithium and of 
some s-process elements are clearly enhanced. Later, 
iKaminskv k. Pavlenko I l|2005f) provided fits to echelle data 
and obtained a similar result, [Fe/H] = —0.3 ± 0.2. These 
results were obtained using a static LTE model and are 
therefore very dependent on the model atmosphere and 
spectrum synthesis assumptions. 

In this paper we discuss fits to the spectral energy 
distribution of V838 Mon using data obtained in 2002 
November. The observational data used in this paper are 
described in Section 2. Section 3 explains some of the back- 
ground to our work and some details of the procedure 
and input parameters used. The results are discussed in 
Section 5. 



2. Observations 

The spectrum we consider was obtained on 2002 
November 6 with the Kast spectograph on the Cassegrain 
focus of the Shane 3-m telescope at Lick Observatory. The 
resolution of the spectrum is A/AA = 3,000. The observ- 
ing and data reduction processes are the same as those in 
IRushton et all l)2005bh and are not repeated here. 



3. Computation of theoretical spectra 

Theoretical spectral energy distributions (SEDs) were 
computed for mo del atmospheres o f gian ts from the 
NextGen g rid of |Hajjgxhj]dt_et aL with solar 

metallicity ijAnders ^r^revesse*TIl989|l . 

In our computations we used a number of model at- 
mospheres, with T cB = 2000-2200 K and log.g = 0,0.5. 
Computations of the synt hetic spectra wer e carried out us- 
ing the program WITA6 l|Pavlenko Il2000h . assuming LTE 
and hydrostatic equilibrium for a one-dimensional model 
atmosphere without sources and sinks of energy. 

The equations of ionization-dissociation equilibrium 
were solved for media consisting of atoms, ions, and 
molecules. We took into account ~ 100 components 
( Pavlenko l200fi) . The constant s for the e quations of chem- 
ical balance were taken from lTsuiil l(l973h . Molecular line 
data were taken from a variety of sources: 



(i) the TiO line lists of IPlez I <ll998h: 

(ii) CN lines from CDROM 1 8 llKuruczlll 993h : 

(hi) CrH and FeH lines from iBurrows et al. I l)2002j) and 
iDulick et al. 1 1 20031). respectively: 

(iv) lines of H^ 6 l|Barber et ai~ll2006|) : 

(v) absorption by VO, and by a few molecules of (in the 
case of V838 Mon) lesser impor tance, was comp uted 
in the JOLA approximation fseelPavlenko 2fJfJfJ); and 

(vi) atomic line list from VALD 1 Kupka etahl 1999ji . 



The profiles of molecular and atomic lines are de- 
termined using the Voigt function H(a,v). Their nat- 
ural (C2) and van der Waals broadenin g (Cq) parame- 
ters are taken from lKunka et al~l lll999|) or. in their ab- 
sence, computed following llTnsoldl Jl955h. Owing to the 



low temperatures in M-star atmospheres, and the conse- 
quent low electron densities, Stark broadening may be ne- 
glected; on the whole the effects of pressure broadening 
dominate. Computations for synthetic spectra were car- 
ried out with a step of 0.003 /im, for microturbulent ve- 
locity Vt = 5kms _1 . 

There are several reasons to suspect that dust might 
be pres ent in the atmosph e re of V 838 Mon, see discus - 
sion in iBaneriee fc Ashokl l(2002l) : lEvans et all l|2003h : 
IRushton et al. I l |2005aft . However, examination of the re- 
sults of preliminary runs of our models without dust show 
that the density of condensibles (such as X, Y, Z) at the 
appropriate distance is low, so that precipitation in the 
photosphere is extremely unlikely. Hence, we do not in- 
clude dust in our models. 

In addition, we do not take account of the effects 
on the model atmos phere of the presence of a B3 star 
l|Munari et al. ll2005|) . Over the wavelength range we are 
modelling (6000-10000 A), the SED is dominated by the 
cool star and the contribution of the B3 star is negligible. 
There may be irradiation effects to be taken into account 
if the hot component is part of a V838 Mon binary system, 
but we do not pursue this here. We will, however, include 
the effects of the B3 star to the overall SED after fitting 
the red end of the spectrum (see Section ^31 below). 

The relative importance of the different opacities con- 
tributing to our synthetic spectra is shown in Fig. ^ 



4. Results 

4.1. Dependence of fits on input parameters 

Changes in our input parameters, i.e. T c g, log g, [Fe/H], 
v t , affect the theoretical spectral distributions (see Fig.[3J) 
in different ways, as follows: 

(i) The response of the spectrum to changes in log g is 
relatively weak due to the absence of strong atomic 
lines with extended wings. Due to low densities in the 
atmosphere of V838 Mon pressure broadening cannot 
be significant. Saturated TiO and VO band heads show 
rather low sensitivity on log g. 

(ii) Due to the same reason, changes in the microturbulent 
velocity mainly affect the heads of strong molecular 
bands. 

(iii) On the other hand, changes in metallicity reduce the 
molecular absorption over the entire spectral region. 

(iv) Due to the low temperatures in the atmosphere of 
V838 Mon changes in the effective temperature af- 
fect the slope of the computed spectra in the region 
0.6 — 1 jum. This effect can be used to determine the 
effective temperature of V838 Mon (see Section E^J). 

Comparison of the computed spectrum with that of 
V838 Mon in 2002 November suggests that the observed 
fluxes at A >0.9 and A <0.8 /jm are better fitted with a 
slightly metal-deficient atmosphere ([Fe/H = —0.5). This 
is in agreement with the analysis by Kipper et al. (2004) 
and Kaminsky & Pavlenko (2004). 
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Fig. 1. Top: contribution of different molecu les to the formation of the spectrum of V838 Mon. The two lower spectra 
represent observed fluxes of L-dwarf Kelul l|Leggett et al. ll2002l) and V838 Mon. Bottom: Identification of the main 
molecular features in the spectrum of V838 Mon. Data for the figure plotting and color version of the plots are available 



on ftp://ftp.mao.kiev.ua/ %2f/pub/users/yp/Fig.l 
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Fig. 2. Dependence of computed spectra on the gravity (top left), microturbulent velocity (top right), metallicity 
(bottom left), and effective temperature (bottom right). The dereddened observed spectrum of V838 Mon is shown by 
the boldfaced line. Model atmospheres are labeled on the figure as T ff/log/[M/H], i.e 2000/0. 0/-1 means T e g= 2000, 
log g =0.0, [M/H] = -1. 



4.2. Comparison with the spectrum of Kelul (L2) 

We compare (see Fig. ^) the observed spectrum 
of V838 Mon with that of the brown dwarf Kelul 
( Legsrett et al. ll2002|l . The L2 brown dwarf Kelul has ap- 
proximately the same effective temperature as V838 Mon 
but, by definition, its log g is very much higher. There are 
several significant points to note: 

(i) Heads of VO bands are clearly seen on at 0.74, 0.79, 
0.87 and 0.96 ^m. 

(ii) The molecular bands of TiO and VO in the spectrum 
of V838 Mon are very intense. Due to the low densities 
these species are practically undepleted in the atmo- 
sphere of V838 Mon despite the low temperatures. On 
the other hand, Ti and V containing species are com- 
pletely depleted in the atmospheres of L-dwarfs. 

(iii) The pressure broadened K I and Na I lines are promi- 
nent in the optical spectra of L dwarfs IPavlenkol 
(2000). Instead of L-dwarfs spectra the alkali lines are 
absent or very weak in the spectrum of V838 Mon. 

(iv) The electronic bands of CrH and FeH are much weaker 
in the spectrum of V838 Mon than they are in Kelul. 



4.3. The effective temperature of V838 Mon in 2002 

A comparison of the observed spectrum of V838 Mon with 
theoretical spectra computed using NextGen model atmo- 
spheres is presented in Fig. |21 in which we use a logarith- 
mic scale for the purpose of illustration as the observed 
fluxes cover several orders of magnitude. The instrumen- 
tal broadening was modelled by Gaussian profiles set to 
the resolution of the observed spectra (FHWM = 2 A). 
The spectrum was d ereddened by E(B — V) = 0.87 mag 
l|Munari et al. Il2005|) . 

The following are worth noting: 

(i) We can fit the overall slope of the spectrum in the near- 
IR/optical region, together with most of the spectral 
features. The spectrum is dominated by TiO bands, 
although some bands of VO are also prominent. 

(ii) This provides clear evidence that, on 2002 November 
6, C/O < 1 in V838 Mon. 

(iii) The dependence on logg is rather weak, as expected 
(see Section 14.10 . Strong atomic lines are practically 
absent in the spectrum. 
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Fig. 3. Fit of theoretical fluxes to the observed spectrum of V838 Mon. Observed and dereddened spectra are shown by 
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Fig. 4. Left: Fit of the theoretical spectrum with T c g = 19000 K and logg = 4.0 to the hot companion of V838 Mon. 
Right: The Balmer jump region of the hot companion spectrum. To simplify the picture we exlude the metal line 
absorption from the accounted opacity sources. 



(iv) The dependence on T e g is much stronger. From com- 
parison of the computed and observed SEDs, we con- 
clude that the effective temperature of V838 Mon was 
T ofr = 2000 ± 100 K in 2002 November. 



Our value of T c ff for 2002 No vember 6 comp ares 
favourably with that determined by iTvlendal (|2005h for 
2002 day 301 (October 28, the nearest in Tylenda's com- 
pilation); he obtained T c ff = 2180 K (formally this corre- 
sponds to a spectral class M8.5). 
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4.4. The hot component in the V838 Mon spectrum 

Our observed spectrum of V838 Mon shows a rise in the 
flux toward the blue end of the spectrum, at A < 6500 A, 
which can be attributed to the presence of the "hot" star; 
this may be a background or foreground object, or even a 
member of a V838 Mon binary system (see discussion in 
iMunari et al. 12005^ . It is worth noting that all the Balmer 
hydrogen lines, from H/3 to the Balmer jump, are clearly 
seen in the observed spectrum. 

We first model the contribution of the hot component 
by a simple Rayleigh- Jeans law F\ = const/ A 4 and nor- 
malise both fluxes to the same value at 0.67 /im. This com- 
paratively simple procedure allows us to fit to the observed 
spectrum over a wide spectral range, from the blue to the 
near-infrared. 

We next computed a set of theoretical spectra using 
the standard Kurucz model and WITA6, for different val- 
ues of T e ff and logg. A comparison of the observed spec- 
trum of the hot companion with theoretical spectra is pre- 
sented in Fig. 0] Unfortunately, our low-resolution spec- 
trum, with its low signal-to-noise ratio in the blue, cannot 
provide a reliable estimate of the atmospheric parameters 
(effective temperature, gravity and/or metallicity) of the 
hot component. From the fit to the observed spectrum we 
estimate T efr = 18000 - 23000 K an d log. q = 4.0. 

From evolutionary calculation ijSchaller et al. I fl992) 
we find the mass of the B star; we then computed the pos- 
sible range of spectral parallaxes for our hot component. 
We obtain M = 6 M , d = 7.2 kpc for T cS = 18000 and 
M = 9 M , d = 14.4 kpc for T off = 23000 K, respectively, 
i.e only in the case of T e s < 21000 K the hot component 
can be a member of the V838 Mon binary system. 



5. Discussion and conclusions 

It seems that the pseudophotosphere of V838 Mon was 
evolving sufficiently slowly that its optical spectrum 
could be approximated b y a sequence of classical, i.e . 
static model atmospheres ( Kaminskv fc Pavlenko ifeoOSh : 
we find that the standard computational procedures for 
SEDs, applied to V838 Mon, give results that are reli- 
able enough to determine the effective temperature. It is 
worth noting that the presence of P Cyg profiles indicates 
expansion velocities up to 200kms _1 . However the emis- 
sion components of the P Cyg profiles form at the outer 
boundary of the envelope, well above the region where the 
observed continuum and absorption lines form. 

There are distinct differences between the SED of 
V838 Mon and of late-type dwarfs with similar spectral 
type, which obviously arise from their differing physi- 
cal nature and the input physical data. In the case of 
V838 Mon we have modelled its SED at an advanced stage 
of its evolution and find that its optical spectrum is formed 
by absorption of the saturated bands of VO and TiO. 

We note that our models do not satisfactorily repro- 
duce the absorption band located around 0.75 /zm; this 
absorption cannot be identified as TiO or VO bands (see 



Fig. Q). Hypothetically, the strong band of CH4 could be 
formed in a low-temperature regime at these wavelengths. 
This band, as well as the bands at 0.619 and 0.890 /im, are 
known to be present in the spectra of Jupiter and other 
gas giants. Possibly the band may be formed in the out- 
ermost layers of the extended envelope of V838 Mon. We 
plan to provide a more detailed study of the molecular 
features in a forthcoming paper. 

Nonetheless, we have obtained a qualitatively good 
fit of the spectrum of V838 Mon over a wide spectral 
range (0.4-1 fim) with a model atmosphere having T e g = 
2000 ± 100 K and C/O < 1, and a hot component consis- 
tent with Teff=18000 - 23000 K and log g =4.0 ± 5; this 
agrees well with previous studies l|Tvlenda I l2005h . The 
cool component is slightly better reproduced with a mod- 
erately sub-solar metallicity model. 
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